Source of Acquisition 
NASA Marshall Space Flight Center 


A ROTATING PLUG MODEL OF FRICTION STIR WELDING 

HEAT TRANSFER 

J. K. Raghulapadu 1 , J. Peddieson 1 *, G. R. Buchanan 2 

Departments of ’Mechanical Engineering and 2 Civil and Environmental Engineering 
Tennessee Technological University 
Cookeville, Tennessee 38505 


A. C. Nunes 

Materials Processes and Manufacturing Department 
NASA Marshall Space Flight Center 
Huntsville, Alabama 35812 


Abstract 

A simplified rotating plug model is employed to study the heat transfer 
phenomena associated with the friction stir welding process. An approximate analytical 
solution is obtained based on this idealized model and used both to demonstrate the 
qualitative influence of process parameters on predictions and to estimate temperatures 
produced in typical friction stir welding situations. 
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Introduction 


The primary purpose of this paper is to develop an approximate closed form 
solution applicable to the description of the heat transfer phenomena associated with the 
friction stir welding (FSW) process. This solution extends classical quasi-steady 
conduction heat transfer solutions for moving sources and hot spots to the modeling of 
FSW heat transfer. A secondary purpose is to use this solution to discuss certain 
qualitative aspects of FSW heat transfer. It is not the purpose of this paper to present 
detail ed pro cess. simulations, . Aninteresting finding .is- that-a well known elementary 
conduction heat transfer solution can be employed to estimate the maximum temperatures 
occurring in FSW. 

A schematic of a typical FSW setup is shown in Fig. 1. Plates A and B are to be 
welded along their common boundary (the weld seam). Plates A and B taken together 
will be referred to as the workpiece. The workpiece has a length L and a width W. 
Typically the workpiece is placed on a backup plate (anvil) and clamped rigidly to 
prevent lateral movement. Figure 1 also introduces a cartesian coordinate system (x,y,z) 
and a cylindrical polar coordinate system (r,_,z) which will both be used in later analysis. 

The weld is performed by a pintool. It traverses the weld seam with, a speed V 
(starting at i and ending at f) while rotating with an angular velocity Q . Figure 2 shows a 
schematic of a typical pintool, which consists of a rotating shaft and a rotating nib as 
depicted. The bottom surface of the shaft, which is in contact with the workpiece, is 
called the shoulder. The workpiece has a thickness H. The nib extends a depth D 
(D H) into the workpiece. If D=H the configuration is called “full penetration” while if 
D<H the configuration is called “partial penetration”. The shoulder radius is R s and the 
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nib radius is R (R<R S ). To begin the process, the nib is brought into contact with the 

weld seam in the vicinity of i (see Fig.l). The contact between the nib and the top 
surface of the workpiece results in the frictional generation of heat. This softens the 
nearby metal and allows the nib to be forced into the workpiece to achieve the 
configuration depicted in Fig. 2 with full shoulder contact. A region of yielded (but not 
melted) metal forms in the immediate vicinity of the pintool which reaches a maximum 
temperature around 80% to 90% of the melting temperature of the workpiece. The 
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pintool moves down the seam the metal cools leaving a completed weld. In the vicinity 
of f (see Fig. 1), the nib runs off the plate (usually onto a runoff tab), is withdrawn 
leaving a hole, or, in the case of the retractable pintool invented at Marshall Space Flight 
Center, is retracted into the shaft without leaving a hole at the end of the seam. 


Some representative papers dealing with thermal and mechanical modeling of 
FSW are those by Chao and Qi [1], Ulysse [2], Hyoe et al. [3], Chao et al. [4], Colegrove 
and Shercliff [5], Song and Kovacevic [6,7], Schmidt et al. [8], Chen and Kovacevic [9], 
Soundararajan et al. [10], Schmidt and Hattel [11], Colegrove and Shercliff [12], Cho et 
al. [13], and Zhang et al. [14]. Several other pertinent publications are referred to therein. 
An inspection of the above cited references shows that there are two basic approaches to 
FSW heat transfer simulations. A considerable portion of the heating which maintains a 
FSW process is provided by mechanical dissipation heating in the yielded region. Thus, 
the most accurate approach is to perform a combined thermomechanical (FSW-TM) 
simulation (see, for instance, [2], [5], and [9-14]). This involves numerical methodologies 
that require extremely fine grids in the pintool vicinity in order to capture the details of 



the small yielded region and, thus, are computationally intensive. Because of this, many 
investigators (see, for instance, [1], [3], [4], and [6-8]) have performed thermal 
simulations (FSW-T) in which the heating is regarded as prescribed. Since the heating 
depends on the details of the mechanical behavior in the yielded region, FSW-T models 
are not as complete as FSW-TM models. Nevertheless, they are desirable for certain 
purposes (preliminary parametric studies, for example) because of their reduced 
computational requirements. The present contribution presents an analytical FSW-T 
mod el which is complementary to the numerical FSW-T models .reported i n the a bove _ 
cited references, It employs a simplified “rotating plug” model of FSW to develop an 
approximate closed form solution describing FSW heat transfer. This solution is used to 
evaluate qualitatively the influence of process parameters on predictions. 

Governing Equations 

The details of a typical actual FSW setup are shown in Figs. 1 and 2. The rotating 
plug FSW heat transfer model is created by introducing several idealizations. These are 
enumerated below. 

First, the pintool is idealized as a nib with D = H (rotating plug). The frictional 
heating associated with the rotating plug surface is adjusted to be representative of the 
frictional heating associated with all surfaces of the actual pintool (shoulder bottom, nib 
side, and nib bottom). The procedure used to accomplish this is discussed in the 
“Numerical Examples” section below. 

Second, the process is treated as steady from the viewpoint of the pintool. The 
justification for this is as follows. Typically, the inequalities R S «W and R S «L both 
hold. Thus, except in the immediate vicinities of i and f, the changes of distance between 
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the pintool and the ends of the weld seam are of secondary importance. Also, classical 
solutions for heat transfer associated with moving hot spots and sources (see, for 
instance, Carslaw and Jaeger [15]) indicate that, after a brief transient period, a steady 
state is established from the viewpoint of the hotspot or source. The same physics is 
assumed to apply in the present situation. Thus the process is viewed from the pintool as 
one in which the pintool rotates with angular velocity U relative to the workpiece which 
moves past it with speed V and the coordinate systems shown in Fig.l are regarded as 
hnving their.jorigins. at. the intersection of t he fixed pintool centerline and the bottom of 
the workpiece. 

Third, because of the inequalities stated in the previous paragraph, it is believed 
that the exact geometry of the workpiece outer boundary in the x,y plane will have a 
minor effect on heat transfer in the vicinity of the pintool. Therefore, in order to facilitate 
closed form solutions, the workpiece is idealized as being infinite in the x,y plane. 

Fourth, the mechanical behavior exhibited by the process is highly idealized. In 
reality the workpiece will exhibit a yielded region close to the pintool and an unyielded 
region away from the pintool. Mechanical dissipation in the yielded region produces 
much of the heating responsible for FSW. In the present work, the yielded region is 
idealized as being coincident with the rotating plug surface and the heating due to 
mechanical dissipation is represented by a heat flux at the rotating plug surface. In 
addition the workpiece motion is idealized as a uniform stream parallel to the weld seam, 
thus neglecting the motion in the yielded region necessary for the metal to flow around 
the plug. 
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The simplifications discussed above produce the rotating plug FSW heat transfer 
model to be employed herein. For this model, the standard heat conduction equation 
(see, for instance, [15]) becomes 

d r 2 T + 5 r T/r + dlT/r 2 + d 2 J = -V( cos(0)3 r T - sin(0)3 e T / r) / a ^ 

when written in cylindrical polar coordinates. The left hand side of Eq. (1) represents heat 
conduction while the right hand side represents heat convection by a uniform stream. In 
Eq. (1), T is the absolute temperature and a is the workpiece thermal diffusivityv The 
corresponding radial, transverse, and axial heat flux components can be written 

q r = -kd r T, q.=-k^T/r, q 2 = -kd z T (2) 

where k is the workpiece thermal conductivity. 

To obtain the heat flux which represents the interaction between the pintool and 
the workpiece, it is assumed that the entire pintool surface is acted on by a transverse 
shear stress at the yield value 6 0 . In the terminology of Schmidt et al. [8], this is a “stick” 

model as opposed to a “slide” model. Experimental data are presented in. [8] suggesting 
indirectly that the “stick” model is superior. 

The mechanical power generated by the moment associated with this yield stress 
acting on an area element dA of the rotating plug surface is 

dp = o 0 RQdA (3) 

which can be converted into an equivalent pintool surface radial heat flux 

q.(R,0,z) = dp/dA = o 0 RQ (4) 
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Combining Eq. (2a) (here and subsequently “a” denoted the first of a set of equations, “b” 
the second, etc.) evaluated at r = R with Eq. (4) produces the expression 

3 r T(R,6,z) = -6 0 RQ/k (5) 

The yield stress is a function of both temperature and strain rate. In the present 
model the pintool angular velocity can be used to characterize the strain rate so it could 
equally well be stated that the yield stress depends on both temperature and angular 
velocity. Since the angular velocity is constant during a FSW operation, the dependence 
on temperature will be focused On herein. Typical maximum FSW temperatures are 
80% to 90% of the melting temperature T M , as mentioned earlier. In order to facilitate 
closed form solutions, it was decided to use the straight line idealization 

6 0 = 6 0>0 (1-T/T M ) (6) 

In the high temperature FSW range. Equation (6) captures both the temperature softening 
effect and the vanishing of the yield stress at the melting temperature in the context of a 
linear model. Results presented by Colegrove and Shercliff [12] suggest that the inclusion 
of temperature softening or some equivalent mechanism is necessary to avoid 
unrealistically high temperature predictions. The quantity 6 0>0 is a parameter having the 

dimensions of stress which would be chosen to make Eq. (6) provide the best 
representation of a given set of yield stress/temperature data in the FSW range at a given 
strain rate. Combining Eqs. (5) and (6) yields 

d r T(R, 0,z) = - 6 0i0 RQ (1 - T(R, 0,z)/T M )/k 
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Because heating is also provided by sliding friction over some portion of the shoulder 
surface (not included explicitly in the present model), some additional adjustment of the 
parameter o 00 will normally be needed in order to obtain realistic predictions. The 

heating model employed herein has many features in common with those discussed in [6- 

8 ], 

It should be mentioned that the closed form solution to be obtained subsequently 
could be readily imbedded in an iteration scheme which would deal with an arbitrary 
yield stress/temperature relation by treating both o 0 0 and T M appearing in Eq. (7) as 

adjustable parameters and recomputing them at each iteration to get the best straight line 
fit to the yield stress/temperature behavior in the vicinity of the temperature calculated at 
the previous iteration. This subject is not pursued further herein. 

The workpiece/anvil and workpiece/top environment interactions will be 
represented herein by the respective heat transfer coefficients h 0 and h H . Then the 
definition of a heat transfer coefficient, combined with the evaluation of the appropriate 


element of Eq. (2) at the appropriate surface leads to 

d z T(r, 0,0) = h 0 (T(r, 6,0) - T A )/k (8) 

0 z T(r,0,H) = h H (T„ - T(r,0,H))/k (9) 

where X*, is the environmental temperature and T A is the anvil temperature. Also 
T(r,0,oo)=T„ (10) 

Equations (1) and (7) - (10) form the boundary value problem 
d 2 J + 3 r T/r + d 2 T/r 2 + d 2 J + P(cos(0)c7 r T T - sin(0)a e T/r)/R = 0 (1 1) 

R 



d r T(R, 0, z) - T(R, 0, z) /(eR) = — T M /(eR) 
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T(r ) e,oo) = T <0 (13) 

d 2 T(r,0,O) -B o T(r,0,Q)/H = -B 0 T A /H (14) 

d z T(r,0,H) + B H T(r,0,H)/H - B H T„ /H (15) 

InEqs. (11)-(15) 

P = VR/a, B H = h H H/k, B 0 = h 0 H/k, e = kT M /(fo 00 QR 2 ) (16) 


are the respective Peclet number, Biot numbers associated with the top and bottom 
workpiece surfaces, and a quantity having no accepted name which will be called the 
“mechanical dissipation number” herein. A correction factor f has been inserted into the 
mechanical dissipation number to account for the fact that the surface area of the rotating 
plug is less than the surface area of the actual pintool. 

It should be noted that the Peclet number is a dimensionless measure of the 
welding translation speed, the Biot numbers are dimensionless measures of the degree of 
thermal interaction between the workpiece and its surroundings, and the mechanical 
dissipation number is a composite dimensionless quantity involving the welding rotation 
speed and the yielding behavior of the workpiece. Thus, the effect of process parameters 
on predictions can be characterized efficiently in terms of these dimensionless numbers. 

Approximate Closed Form Solution for Small Peclet and Biot Numbers 

Inspection of the various parametric combinations discussed in the previously 
cited literature reveals Peclet and Biot numbers in the respective ranges 0. 1-0.3 and 
0.001-0.02. This information suggests that approximations based on the smallness of the 
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Peclet and Biot numbers should be useful. Further, the range of Biot numbers suggests 
equating the Biot numbers to zero in Eqs. (14) and (15). This creates boundary 
conditions describing insulated top and bottom workpiece surfaces and allows a solution 
of the form T = T(r,e) which reduces Eqs . (1 1), (12), and (13) to the respective forms 

d 2 r T + d r TVr + 5 2 T/r 2 = -P(cos(e)d r T - sin(e)3 J/r)/R 

3 r T(R,e) - T(R,e:)/(aR) - -T M /(SR), T(°°,e) = T M (17) 

Equations (17) describe a two dimensional version of the FSW-T problem. Examples of 
other recent two dimensional FSW simulations are provided by the papers of Cho et al. 
[13] and Zhang et al. [14]. 

A simple approximate closed form solution of the problem defined by Eqs. (17) 
can be found for small Peclet numbers using an approach similar to that discussed by 
Trivedi and Srinivasan [16]. Omitting the details for the sake of brevity and using the 
small argument form of a Bessel function, the result can be written in the reduced 
temperature form 

(T-TJ/(T m - TJ = ip (r) exp(- Pr cos(0) / (2R)) (18) 

where 

rp = K 0 (Pr/2R))/(e + ln(2/P)) (19) 

In Eq. (19) Ko is a zeroth order modified Bessel function of the second kind. Near the 
pintool the small argument form of _ can be used to further simplify Eq. (18) to 
(T-T.)/(T M -T.)-^ 0 (r)exp(-Prcos(e)/(2R)) - (20) 

where 



t|) 0 = ln(2R/(Pr))/(s + ln(2/P)) 


( 21 ) 


The combination of Eqs. (20) and (21) with the exponential term omitted is just 
the reduced temperature associated with axisymmetric steady heat conduction in a 
cylindrical shell having respective inner and outer radii 

r i - Rj r o = 2/P (22) 

and respective inner and outer temperatures 

Tj = T„ + ,(T M - T„ )ln(2 / P) /(e + ln(2 / P)), T 0 = T M (23) 

(see, for instance, Incropera and DeWitt [17]). Thus, this elementary solution can be 
adapted to the description of the temperature distribution near the pintool in FSW. 

Evaluating Eq. (21) at the pintool surface yields 
<j>-l/(l+t) (24) 

where 

<() = ^ 0 (R), t, = e/ln(2/P) (25) 

It can be shown that, to the order of approximation inherent in Eq. (18), <p represents the 
average reduced temperature on the pintool surface. For 0-.l<P<0.3 the maximum and 
minimum reduced temperatures will differ from the average by a maximum of about 
15%. The elementary cylindrical shell solution can, therefore be used to estimate the 
maximum temperatures produced in a FSW operation. 

The combination of Eqs. (18) and (19) can be rewritten as 

CT - T„)/(T m -T„) = QK 0 (Pr/(2R))exp(-Prcos(8)/(2R))/(2jck) (26) 
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This is identical to the classical solution for quasi-steady conduction due to a moving line 
source of heat (see, for instance, Carslaw and Jaeger [15]) with source strength (power 
per unit length) 

Q = 2jtk(T M - ) /($ + ln(2 / P)) (27) 

Thus, except in the immediate vicinity of the pintool, the workpiece cannot distinguish 
between mechanical dissipation and line sources of heat. 

The combination of Eqs. (18) and (19) can also be written 

"(f-i l)/(i^-X>K^ ( 28 ) 

For large distances from the pintool, the large argument limit of Ko can be used to 
simplify Eq. (28) to 

(T-TJ/(T m -T m ) = (jcR/(P(x 2 +y 2 ) 1 / 2 )) 1 / 2 exp(-P((x 2 +y 2 ) 1/2 
+ x)/(2R))/(e + ln(2/P)) 

Equation (29) indicates exponential decay of the recuced temperature along lines parallel 
and perpendicular to the weld seam except along the weld seam itself behind the pintool 
where the decay is algebraic like inverse square root of x. This is a very slow decay and 
indicates the need for very long computational regions behind the pintool when numerical 
FSW-T simulations are performed. 

All results given thus far represent the viewpoint of an observer moving with the 
pintool relative to a fixed workpiece (or, equivalently, describe a workpiece moving 
relative to a fixed pintool). The viewpoint of a fixed observer at a distance X to the right 
of point i (see Fig. 1) on a fixed workpiece can be created by replacing x in Eq. (28) by 
X-Vt (t being the time since process initiation) to get 
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(T-TJ/(T m -TJ = K 0 (P((X- Vt) 2 + y 2 ) 172 /(2R))exp(-P(X 
-Vt)/(2R))/(s + ln(2/P)) 

This describes a rapid increase in temperature as the pintool approaches followed by a 
slower decay as the pintool departs. This is in qualitative agreement with the 
experimental data presented in several of the previously cited references (see, for 
instance, figure 9 of [1], figure 13 of [7], figure 4 of [9], and figure 10 of [1 1]). 

The simple form of Eq. (24) makes possible the definition of a sensitivity factor 

„S.=.^d(j> /-<j>) l(d% I £) = (£/cj>)(d<j>./dQ = £/<!+.£) - (31) 

This can roughly be thought of as the ratio of the infinitesimal percentage change in (j) to 
the infinitesimal percentage change in _ producing it. The sensitivity factor S is an 
increasing function of _ with a minimum value of zero and a maximum value of unity. A 
small percentage change in _ will, therefore, always produce a smaller percentage change 
in (j). Thus, (j> is not highly sensitive to small changes in _. This allows some freedom in 
choosing the fitting parameter _. As a specific example, the parametric configurations to 
be discussed in the next section suggest that _ seldom exceeds 1/2. The value _ =1/2 
corresponds to S=l/3, indicating that a small percentage change in _ in the vicinity of 1/2 
will produce a percentage change in <J> one third as large. 


Numerical Examples 

Since the forms of the results presented in the previous section are so simple, it 
seem unnecessary to report extensive tabular or graphical representations of predictions. 
Nevertheless, a few numerical examples will be presented in this section with the purpose 
of investigating the issue of tuning. As mentioned earlier, no FSW-T model is fully 



determinate because the heating cannot be completely characterized without a knowledge 
of the mechanical behavior. For this reason, tuning is a feature of all FSW-T models (see, 
for instance, [1-4] and [6-8]). Several choices of tuning parameters have appeared in the 
literature. In the present simplified model there is no choice. The quantity _ 0>0 must be 
selected. It is of interest, therefore, to determine the sensitivity of predictions to _ 0 , 0 and 
to see if a value chosen on the basis of one set of experimental data can predict other sets 
of experimental data. 

To deal with realistic configurations-, it is necessary to select the correction factor 
f to account for the actual pintool geometry shown in Fig. 2. There is no unique way to do 
this. Here it is accomplished by assuming that a constant transverse shear stress 'at the 
yield value is acting on the bottom and side surfaces of the nib and the bottom surface of 
the shoulder and finding the equivalent heat flux associated with the mechanical power 
due to the corresponding moment about the pintool centerline (similar to the 
methodology employed in [8]). Doing this shows that Eq. (4) must be replaced by 

q n =6 0 rQ (32) 

to take into account the variable moment arms on the nib and shoulder bottoms. In Eq. 
(32) q n is the heat flux normal to the pintool surface. According to this model, the net 
heating rate due to the mechanical dissipation is 

Q - ^q,dA (33) 

where A is the combined nib bottom and side areas and the shoulder area. Carrying out 
this integration and defining an equivalent radial heat flux for the rotating plug model as 
q r (R,e,z) = Q/(25RII) (34) 
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(the denominator of Eq. (34) being the rotating plug surface area) produces Eq. (4) with 
the o 0 appearing therein increased by a multiplicative factor of 

f =D/H + (R/H)(R s /R) 3 /3 (35) 

For the full penetration configuration (D — H) the corresponding result is 
f - 1 + (R / H)((R S / R) 3 — 1)/3 (36) 

Four typical configurations based on the previously cited literature are 
characterized in columns 1-7 of Table 1. The values of k and a correspond to 6061 A1 
which is the workpiece material in all five of these configurations. The value of o 0 0 = 

60.0 MPa — 60.0 N /mrn 2 was obtained by timing to make the model predict temperatures 
in the FSW range. The need for such tuning in all FSW-T models was discussed earlier. 
Obviously, tuning can be done in many ways and no attempt to optimize the process was 
made herein. 

Columns 8-11 of Table 1 show the values of pertinent intermediate quantities and 
column 12 shows the average reduced temperatures (<j>) at the pintool surface. Column 13 
contains the corresponding values of the predicted pintool surface average absolute 
temperature (T). The values reported in column 1 1 provide support for the estimate _<l/2 
used in the previous section. 

Column 14 of Table 1 presents estimates of the observed maximum absolute 
temperatures (T) in the experiments discussed in each of the references cited. Estimates 
are required both because it is possible to measure temperatures in only a limited number 
of locations in a FSW experiment (and practically impossible to measure the maximum 
temperature) and because not all data are reported in publications. Even allowing for 



some estimation error, it can be said that the agreement between predictions and 
observations is remarkably good, considering the simplicity of the model. It should be 
emphasized that a single value of _ 0)0 has been used to make all five predictions. No 
attempt has been made to correct for the dependence of yield stress on strain rate or the 
fact that part of the pintool surface may be subjected to sliding friction. 

Conclusion 

.In .the , forgoing, a simplified rotating plug, model of friction stir wel ding heat 
transfer was used to develop a relatively simple closed form approximate solution 
predicting the workpiece temperature distribution produced by FSW. It was found that 
the temperature distribution in the immediate vicinity of the pintool could be 
characterized by a well known elementary conduction heat transfer solution. It was 
further observed that the average reduced temperature at the pintool surface depended on 
only one composite dimensionless parameter containing information about the 
translational and rotational weld speeds, the pintool geometry, and the workpiece 
material properties. The sensitivity to this parameter was found to be moderate. A single 
value of a tuning parameter was found to be able to produce predictions in good 
agreement with data from four different sets of FSW experiments involving 6061 A1 as 
the workpiece material. This suggests that once the model is tuned using one set of data 
for a particular workpiece material, it should be able to deal with other configurations 
pertinent to the same workpiece material. It is believed, therefore, that the present 
rotating plug FSW-T model should be useful for quick preliminary parametric studies of 
proposed FSW configurations. 



Acknowledgement 

This research was supported by the National Aeronautics and Space 
Administration George C. Marshall Space Flight Center through Grant/Cooperative 
Number NCC8-223 for the National Center for Advanced Manufacturing - Louisiana 
Partnership. Any opinions, findings, and conclusions or recommendations expressed in 
this material are those of the authors and do not necessarily reflect the views of the 
National Aeronautics and Space Administration. 

Nomenclature 

A=p intool area 

B-Biot number 

D=pin depth 

f=area correction factor 

H=workpiece thickness 

h=convection heat transfer coefficient 

Ko=zeroth order modified Bessel function of second kind 

L=workpiece length 

k=thermal conductivity • 

P=Peclet number 
p=power 

Q=heat source strength 
q=heat flux 
R=pin radius 
Rs=shoulder radius 
r=polar coordinate 
rpshell inner radius 
r 0 =shell outer radius 
S=sensitivity factor 
T=temperature 
TM^melting temperature 
Ti- shell inner temperature 
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T 0 =shell outer temperature 
T A =anvil temperature 
T ^environmental temperature 
t=time 

V=pintool velocity 
W=workpiece width 
X=cartesian coordinate 
x=cartesian coordinate 
y=cartesian coordinate 
z=cartesian or polar coordinate 

_=thermal diffusivity 
_=mechanical dissipation number 
_=composite dimensionless parameter 
_=polar coordinate 
_o=yield stress 
_o,o=fitting parameter 
<j>~ average reduced pin temperature 
_=reduced temperature function 
_ 0 =shell reduced temperature function 
_=pintool angular velocity 

subscripts 
H=workpiece top 
n=normal 

o=workpiece bottom 
r=radial 
z=axial 
=transverse 
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Figure captions 

Fig. 1 Schematic of typical friction stir welding configuration 

Fig. 2 Schematic of a typical pintool (partial penetration configuration) 








Table 1 Maximum temperatures for typical parametric combinations involving 606 1A1 
(k = 170 N/(sK),a = 61mm 2 /s,T M = 855K 

t 0 0 = 60.0N7 mm 2 , T„ = 300K) 
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